N95- 23313 


Aeroacoustic computation of gust-blade interaction 

by 

James E. Martin 
Department of Mathematics 
Christopher Newport University 
Newport News, VA 23666 


To better understand and address the challenges faced in computing the acoustics of 
flow fields, test problems must be considered. In the present study, the sound radiated by 
the interaction of a flat plate with an oncoming gust containing a two component, mean 
velocity is computed. The gust has a uniform mean flow in x with Mach number A/oo 
equal to 0.5. The gust’s mean velocity in y is of smaller amplitude and is given by 


v = 0.1 sin 
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This problem has been posed for an upcoming ICASE/LaRC workshop on benchmark 
problems in computational aeroacoustics. 

A plate with a length of 30 units in x is used. The plate is assumed to be in- 
finitesimally thin and is centered at the origin. All variables are made dimensionless 
using the scales specified. Acoustic quantities are obtained by numerically integrating 
the linearized Euler equations. Integration is performed on the computational domain 
-100.0 < x < 100.0, -100.0 < y < 100.0, using unit length grid spacing in x and in y. 

An integration scheme is sought which will provide accurate solution to the small quanti- 
ties of interest at a minimal computational expense. Results indicate that with the given 
discretization a scheme of minimal fourth order accuracy might be adequate to approxi- 
mate the waves within the given flow. Thus, a variation of the MacCormack scheme with 
fourth order accuracy in space and second order accuracy in time was chosen. A scheme 
with sixth order accuracy in space has also been implemented and results compared with 
those of the fourth order accurate scheme. 

To ensure no mass flux, zero normal velocity is assigned at the plate. This condi- 
tion will induce a discontinuity in the pressure across the plate location. Values for the 
perturbation pressure p f along the surface of the plate are obtained using a one-sided, 
third order Taylor expansion, such that p f y = 0. In accordance with the Kutta condition, 
perturbation pressure at the trailing edge is assigned to zero. In the far field, radiation 
boundary conditions have been implemented. The effectiveness of the far field conditions 
are validated by computing in a larger computational domain and comparing the results. 

Early in time, sound waves begin to radiate from the plate. A Doppler effect is 
observed. After the initial transients disappear, the strongest waves leave the trailing edge 
at an approximate 45 degree angle. The intensity pattern of pressure fluctuations shows 
five lobes (of increasing magnitude with increasing downstream direction) emerging from 
the plate. Undesirable short wave contaminants are observed in the computed pressure 
distribution along the plate surface. For a more accurate solution at small scales, a more 
refined discretization will be required. 
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